Electrochemical production of ammonia and catalyst therefor

ABSTRACT

An iron-containing Chevrel phase material, contains iron and Mo 6 S 8  clusters, in particular an iron-containing Chevrel phase material having a formula Fe x Mo 6 S 8 , wherein 2≤x≤4. The iron-containing Chevrel phase provides an efficient catalyst for the electrochemical production of ammonia from water and nitrogen gas.

BACKGROUND

Ammonia (NH₃) is an essential feedstock for many industrial processes including agriculture, chemical production and pharmaceuticals.¹ It is also one of the most promising carriers for renewable electricity as it has a high energy density of 3 kWh kg⁻¹ and a 17.7 wt % hydrogen content.² The demand for ammonia has been continuously increasing and it is currently among the most highly produced inorganic chemicals (at ˜500 million tons per year).³ Unfortunately, the industrial production still relies on the Haber-Bosch process that was invented in the 1900s. This process requires centralized plants and operates under harsh conditions using hydrogen derived from fossil fuels (at 300˜400° C., ˜250 atm), it consumes 1˜2% of the global energy supply and is responsible for >1% of total greenhouse gases emissions.^(4, 5) Therefore, it is crucial to develop alternative technologies for more affordable and sustainable NH₃ production, and the electrochemically based catalytic systems are among the most attractive candidates.⁶⁻⁸ Such systems could in principle operate under ambient conditions and are integrable with renewable electricity to directly produce ammonia from humidified air in a carbon-neutral manner without reliance on fossil fuels.^(9, 10) It could also be implemented as stackable modules for on-demand and decentralized NH₃ production, therefore substantially mitigating the urgent energy and environmental challenges. ^(5, 11)

The ambient electrocatalytic activation and conversion of N₂ in aqueous electrolytes faces substantial challenges, particularly in the lack of suitable electrocatalysts that are tuned specifically for highly selective NH₃ production.^(9, 12, 13) The past decades of research examined a plethora of nanostructured catalysts, including noble and transition metals such as Ru, Pd, Au, Fe and Ni, doped porous carbon and a variety of metal oxides, nitrides and sulfides.^(12, 14-20) Unfortunately, almost all of these catalysts are plagued by slow kinetics and low Faradaic efficiency, with NH₃ production rates that are far from being competitive with the Haber-Bosch process.²¹ Furthermore, the conversion to NH₃ involves transferring 6 electrons and protons to N₂ via complicated, multi-step processes which are very difficult to modulate due to inherent limitations of existing electrocatalysts. This is because the active sites of these catalysts generally have similar binding characteristics and could be best described as single site catalysts;^(14, 15, 22-25) and during N₂ conversion these sites were most likely to first absorb the H-donating species (H₃O⁺ or H₂O, depending on the pH) under electrochemical potential due to the low concentration and the highly inert N≡N triple bond of N₂. This step is likely followed by reaction with N₂ and generating the key *N₂H intermediate for the final production of NH₃, either via proton-coupled electron transfer or hydrogenation with electrochemically generated H_(ad)* processes.^(15, 26) Unfortunately, both steps have high energy barriers and low probability of occurrence so the majority the H-species proceed with the hydrogen evolution reaction (HER), thereby significantly compromising NH₃ selectivity.²⁷ Although these observations suggest that suppressing the HER is essential, it would be more important to design multi-active site catalysts that have separate binding sites for N₂ absorption, proton activation and simultaneous generation of activated H_(ad)* and N₂*.^(9, 28) The synergy of these binding sites may substantially accelerate NH₃ formation when these intermediates are closely oriented on the catalyst surface. Design of such catalysts is an ongoing challenge and necessitates a better understanding and control of the catalytic materials and the reaction mechanism.

In nature, metalloenzymes such as MoFe nitrogenase convert N₂ to NH₃ under ambient condition. These enzymes typically have a N₂ binding protein (the cofactor) and a reducing protein, both containing several subunits with multiple Fe, Mo and S atoms at different oxidation states.²⁹ These subunits provide binding sites that synergistically assist N₂ adsorption and transformation to NH₃.³⁰ These natural biocatalysts have inspired intensive efforts on designing molecular analogues to mimic the process, but rarely on heterogeneous catalysts that are potentially more durable for practical applications.³¹

SUMMARY

In a first aspect, the present invention is an iron containing Chevrel phase material, containing iron and Mo₆S₈ clusters.

In a second aspect, the present invention is an electrode, comprising (1) a conductive substrate, and (2) an iron-containing Chevrel phase material, on the conductive substrate.

In a third aspect, the present invention is a system for producing ammonia electrochemically from N₂ and water, comprising (a) a working electrode, comprising (1) a conductive substrate, and (2) the iron-containing Chevrel phase material, on the conductive substrate. The system also comprises (b) a counter electrode, and (c) an ion-conductive separator between the working electrode and the counter electrode.

In a fourth aspect, the present invention is a facility for manufacturing ammonia, comprising (A) a system for producing ammonia electrochemically from N₂, and water, (B) a power source, (C) a water feed, and (D) an ammonia dispenser.

In a fifth aspect, the present invention is a method of producing ammonia electrochemically from N₂ and water, with a system for producing ammonia electrochemically from N₂ and water, comprising passing electricity through the working electrode and the counter electrode, to produce ammonia from N₂ and water.

Definitions

Faradaic efficiency, FE, of an electrochemical process for producing ammonia from N₂ and H₂O means the number of electrons consumed in the reaction N₂+6H⁺+6e⁻→2NH₃, divided by the total number of electrons consumed by the system, and may be calculated using the following equation:

FE(NH₃)=[3F×c(NH₃)×V]/Q

where F is the Faraday constant (96485 C mol⁻¹), Q is the total charge passed through the electrode, V is the volume of the electrolyte and c(NH₃) is the quantified ammonia concentration. Alternatively, “c(NH₃)×V” in the equation may be replaced by the total number of moles of ammonia produced.

The mass-normalized yield rate of NH₃ for a catalyst used in an electrochemical process for producing ammonia from N₂ and H₂O may be calculated using the following formula

Yield Rate_(mass)(NH₃)=[17×c(NH₃)×V]/(t×m)

where t is the electrolysis time, m is the loading mass of the catalyst, V is the volume of the electrolyte and c(NH₃) is the quantified ammonia concentration. Alternatively, “c(NH₃)×V” in the equation may be replaced by the total number of moles of ammonia produced.

The term “particle size” of a particle means the diameter of a circle having the same area as that of a particle when viewed by transmission electron microscopy. The term “average particle size” mean the average of the particle sizes of a collection of particles.

As used herein, the term “Chevrel phase material” means a compound that contains Mo₆S₈ clusters, such as those illustrated in FIG. 1. A Chevrel phase material may be identified by a variety of methods, including X-ray crystallography and powder X-ray diffraction (XRD).

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic illustration of the crystal structure of Fe₂Mo₆S₈ showing interconnected Mo₆S₈ clusters and the two possible Fe ions sites around the trigonal axis.

FIG. 2 illustrates a working electrode,

FIG. 3 shows a system for producing ammonia electrochemically.

FIG. 4 shows a facility for manufacturing ammonia.

FIG. 5 illustrates the structural characterizations of the as-synthesized Fe₂Mo₆S₈ electrocatalyst by X-ray diffraction pattern.

FIG. 6, FIG. 7 and FIG. 8 are a scanning electron microscopy (SEM) image and transmission electron microscopy (TEM) images of as-synthesized Fe₂Mo₆S₈ at different magnifications.

FIG. 9 is a graph comparing nitrogen absorption-desorption isotherm of Fe₂Mo₆S₈ and Mo₆S₈. The measured BET specific surface area is also included.

FIG. 10 is graphs showing linear sweep voltammetry curves of the electroreduction of N₂ to NH₃ on the Fe₂Mo₆S₈ catalyst in Ar- and N₂-saturated electrolytes, with a scan rate 5 mV/s.

FIG. 11 illustrates the corresponding Tafel plots for FIG. 10.

FIG. 12 are chronoamperometric curves at different electrode potentials (vs. reversable hydrogen electrode (RHE)) for 2 hours of the electroreduction of N₂ to NH₃ on the Fe₂Mo₆S₈ catalyst.

FIG. 13 is the UV-Vis absorption spectra of NH₄ ⁺ quantification using the indophenol blue method.

FIG. 14 is a calibration curves based on the absorbance at 665 nm for NH₄ ⁺ quantification using the indophenol blue method.

FIG. 15 shows a comparison of Faradaic efficiency (curve) and corresponding NH₃ yield rates (bars) under different applied potential for the electroreduction of N₂ to NH₃ on the Fe₂Mo₆S₈ catalyst.

FIG. 16 is a graph showing the durability of the Fe₂Mo₆S₈ catalysts as revealed from the relationship between ammonia content and electrolysis duration for 50 hours (at 0.2V vs. RHE) during the electroreduction of N₂ to NH₃ on the Fe₂Mo₆S₈ catalyst.

FIG. 17 shows the NMR analysis of electrolyte solutions from ¹⁴N₂ and ¹⁵N₂ feeding gases.

FIG. 18 is a survey XPS spectrum of the Fe₂Mo₆S₈ catalyst. No features corresponding to N₁s (squared-off area) was observed, confirming no N species were present in the catalyst.

FIG. 19 is a comparison of UV-Vis absorption spectra of mixed solutions containing 5 ml N₂H₄ coloring solution with either 5 ml fresh or electrolyzed electrolyte. No characteristic absorption peak was observed at ˜455 nm, confirming no formation of N₂H₄ during nitrogen conversion reaction.

FIG. 20 is a schematic illustration of the charge transfer of Fe to Mo₆S₈.

FIG. 21 is the Raman spectrum of Mo₆S₈.

FIG. 22 is the Raman spectrum of Fe₂Mo₆S₈.

FIG. 23 is the Mössbauer spectroscopy of Fe₂Mo₆S₈.

FIG. 24 and FIG. 25 are high resolution X-ray photoelectron spectra of Mo 3d and S 2p.

FIG. 26 and FIG. 27 illustrate the powder X-ray diffraction patterns of Mn₂Mo₆S₈ and Cu₂Mo₆S₈, respectively, confirming their pure Chevrel phase structure.

FIG. 28 is a comparison of N₂ to NH₃ conversion efficiency for Mo₆S₈ and M₂Mo₆S₈ (M=Fe, Mn and Cu) electrocatalysts.

FIG. 29 are LSV curves for hydrogen evolution of Mo₆S₈ and M₂Mo₆S₈ (M=Fe, Mn and Cu) electrocatalysts in Ar-saturated electrolyte.

FIG. 30 is a Tafel plot of Mo₆S₈ catalysts in Ar-saturated electrolyte.

FIG. 31 shows a solid-state system for producing ammonia electrochemically.

FIG. 32 is the X-ray powder diffraction pattern of Fe₄Mo₆S₈.

FIG. 33 and FIG. 34 are SEM images of Fe₄Mo₆S₈.

FIG. 35 is a plot of the linear sweep voltammetry of Fe₄Mo₆S₈ in Ar and N₂ saturated electrolytes.

FIG. 36 illustrates a comparison of Faradaic efficiencies for ammonia production and average ammonia yield of Fe₄Mo₆S₈ at different applied currents.

FIG. 37 illustrates a comparison of Faradaic efficiencies for ammonia production and average ammonia yield of Fe₄Mo₆S₈ as observed within 10 cycles.

DETAILED DESCRIPTION

The iron-containing Chevrel phase chalcogenides, including Fe₂Mo₆S₈ and Fe₄Mo₆S₈, are highly efficient electrocatalysts for selective electrochemical conversion of N₂ to NH₃. With the preferred Fe₄Mo₆S₈, stable Faradaic efficiencies of up to 25% were observed at −0.20 V vs. RHE together with a high rate of 102 μg h⁻¹ mg_(cat) ⁻¹ for NH₃ production in aqueous electrolyte. The formation of NH₃ from N₂ was confirmed from a series of control experiments including the ¹⁵N₂ isotope labeling test, and the catalyst exhibited outstanding stability for at least 50 hours. The intrinsic activities were attributed to the unique atomic configurations of Fe, Mo and S in the iron-containing Chevrel phase chalcogenides, including Fe₄Mo₆S₈, that provide separate but synergistic binding sites for N₂ and H addition. Experimental evidence points to the Fe/Mo sites being responsible for absorbing and activating N₂ with the Fe-promoted S sites providing stronger S—H binding that effectively suppressed the hydrogen evolution reaction. The distinctive coordination environment in the Fe—Mo₆S₈ framework ensures synergy of these active sites that accelerate association of key intermediates for selective NH₃ production.

The iron-containing Chevrel phase compounds may be produced by mixing the elements and/or compounds containing the elements, in the desired proportions, together with grinding, followed by high-temperature heating under an inert atmosphere.³³ Molybdenum and sulfur are always present in a ratio of 6:8. Iron is always present, in an iron:molybdenum:sulfur ratio of x:6:8 where 0<x≤6, preferably 2≤x≤4, and more preferably x=4. In the case where x=1 or 6, the product may contain impurities. Other metals may also be present, based on the formula Fe_(x)M_(y)Mo₆S₈, wherein M is at least one metal selected from the group consisting of elements of Groups 1-15 and the Lanthanide series, 0<x≤6, 0≤y<6, and x+y=1 to 6, preferably Fe_(x)Mo₆S₈ (where y=0 in the more general formula); more preferably M is at least one metal selected from the group consisting of elements of Groups 1, 2, 3, 7, 10, 11, 12, and the Lanthanide series, 2≤x≤4, and x+y=2 to 4; most preferably M is a metal selected from Cu, Cd, Na, Mn and Zn, 2≤x≤4, and x+y=2 to 4. Examples include FeMo₆S₈, Fe₂Mo₆S₈, Fe₃Mo₆S₈, Fe₄Mo₆S₈, Fe₅Mo₆S₈ and Fe₆Mo₆S₈, as well as Fe_(x)Mo₆S₈ having non-integer value of x in between. Iron-containing Chevrel phase compounds of the formula Fe_(x)Mo₆S₈ or Fe_(x)M_(y)Mo₆S₈, where x<2 or 4; or x+y<2 or 4, respectively, may also be formed by forming the phase Fe_(x)M_(y)Mo₆S₈, where x+y=2, with M including Cu or Na, followed by oxidation to cause deintercalation of the Cu or Na; oxidation may be carried out chemically (for example, with a mixed acid containing 8.0 M HCl and 0.05 M HNO₃) or electrochemically. Similarly, for composition Fe_(x)M_(y)Mo₆S₈, where x+y>2, a compound of the formula Fe_(x)M_(y)Mo₆S₈, where x+y=2 may be formed, followed by intercalation with metal ions (such as Li, Na or Cu) through chemical or electrochemical reduction. Examples of metal M in any of the above formulas include the elements magnesium, aluminum, calcium, scandium, chromium, manganese, nickel, cobalt, copper, zinc, gallium, yttrium, ruthenium, rhodium, palladium, silver, cadmium, indium, tin, lanthanum, cerium, praseodymium, neodymium, samarium, europium, gadolinium, terbium, dysprosium, holmium, erbium, thulium, ytterbium, lutetium, rhenium, osmium, iridium, platinum, gold, mercury, lead and mixtures thereof. In an alternative of any of the above compositions, Fe_(x)Mo₆S₈ where x=1, x=1 or 2, or 1≤x≤2, are excluded from the catalyst compositions.

When used as a catalyst, the iron-containing Chevrel phase material may be used as produced, or may be milled to provide a greater surface area. When used as a catalyst, the iron-containing Chevrel phase material may be formed into a catalyst ink, containing the iron-containing Chevrel material, a liquid (such as water, ethanol, acetone or mixtures thereof), a binder (such as a hydrogen-conductive polymer or other water-insoluble polymer), and optionally a conductive material (such as carbon or graphite).

As illustrated in FIG. 2, a working electrode, 10, includes catalyst, 12, on and in electrical contact with a conductive substrate, 14. The conductive substrate may be formed from any conductive material which will not significantly react with the catalyst or the electrolyte solution under the condition of use. Examples of conductive materials include carbon (such as graphite or glassy carbon), steel, nickel, copper, gold, platinum, and alloys thereof. The working electrode may have any convenient shape, such as a rod, disc, or may form the container which holds the electrolyte solution. The catalyst may be attached to the substrate using an adhesive, for example a hydrogen-ion conductive polymer solution, such as a solution of NAFION™ in water and isopropanol. Typically, the catalyst in particulate form is mixed with a liquid adhesive to from a catalyst ink, which is then applied to the surface of the conductive substrate. The electrode is formed when the adhesive dries. In addition, polyvinylidene difluoride (PVDF) and polytetrafluoroethylene (PTFE) may also be used as binders for holding the catalyst to the conductive substrate. In the case of a solid-state system (described below) the adhesive or binders may be used to hold the catalyst on the solid-state electrolyte and/or the conductive gas diffusion layer.

As illustrated in FIG. 3, a system for producing ammonia electrochemically, 20, includes a cathode compartment, 22, and an anode compartment, 24, in ion-conductive communication with each other (during operation), with each holding a cathode electrolyte solution, 26, and an anode electrolyte solution, 27, respectively. A separator, 28, allows ions to pass between the two compartments, but otherwise prevents bulk mixing; examples of a separator may be a hydrogen-ion conductive membrane such as NAFION™ membrane, or simply a porous material which allows ions to pass between the fluids, such as a porous membrane, for example cellulose, cardboard or a glass frit. Within the cathode compartment are the working electrode, 30, and the optional references electrode, 34, which are in contact with the cathode electrolyte solution and in electrical communication with an electrochemical power supply (which may also include an electrochemical monitor, controller and/or work station), 36. Similarly, the counter electrode, 32, is within the anode compartment and is in contact with the anode electrolyte solution. Also illustrate are optional inlets and outlets for the addition or removal, respectively, of reactants and products, from the cathode and anode compartments including a nitrogen gas inlet, 38, an ammonia outlet or cathode electrolyte solution outlet and/or inlet, 40, a water inlet or anode electrolyte solution inlet and/or outlet, 44, and an oxygen gas outlet, 42.

FIG. 31 shows a solid-state system, 100, for producing ammonia electrochemically. The system includes a cathode conductive gas diffusion layer, 108, and anode conductive gas diffusion layer, 110, in contact with the cathode conductive base, 112 (cathode), and anode conductive base, 114 (anode), respectively. A catalysts layer, 104, is in contact with the cathode gas diffusion layer and the solid electrolyte, 102. Also shown is an optional anode product layer, 106, which is in contact with the solid electrolyte and the anode gas diffusion layer; if the anode product layer is absent, then the solid electrolyte is in contact with the anode gas diffusion layer. The anode product layer may contain an oxidation product catalyst. A magnified view of the catalytic layer, 118, shows an intimate mixture of a catalyst, 122, and an adhesive or binder, 124. The cathode conductive gas diffusion layer together with the catalyst layer forms the working electrode, while the anode conductive gas diffusion layer together with the optional anode product layer forms the counter electrode. As shown, nitrogen gas and water (reactants) are supplied to the working electrode, and ammonia is produced, while water is supplied to the counter electrode where oxygen is produced. A carrier gas, such as nitrogen, argon or air may be included with the reactants and water supplied to the electrodes. Liquid water may also be present in the catalyst layer and/or in the gas diffusion layers. A power supply, 116, is present to provide electrical power to drive the production of ammonia. The materials used and the system as a whole are similar to those of a solid-state water electrolysis cell, except using the catalyst of the present application in the catalytic layer, and using a proton conductive material, such as a NAFION™ membrane, as the solid electrolyte. See, for example, ref. 56.

Ammonia may be from by passing an electric current through the working electrode (cathode) and the counter electrode (anode), producing ammonia from water and nitrogen gas, and producing oxygen at the counter electrode. If an electrolyte solution is present, the ammonia produced will typically dissolve in the electrolyte solution, forming ammonium salts. Preferably, the electrolyte solution is acidic, having a pH of less than 7, more preferably a pH of less than 5, for example a pH of 3-4. As the ammonia is formed, it will cause the electrolyte solution to become less acidic, so inclusion of a buffer to help maintain the desired pH is preferred. Over time, a continuous synthesis of ammonia may be carried out by removing the ammonia (or ammonium salt) containing electrolyte, and adding fresh electrolyte solution. The electrolyte at the counter electrode may become more acidic, so addition of acid, water, and/or fresh electrolyte may be desirable. Ammonia may be isolated as ammonium salts from the electrolyte solution, or may be obtained as ammonia gas using chemical methods. In a solid-state system, continuous production of ammonia is also possible, with a continuous supply of reactants and removal of products.

Nitrogen will be consumed at the working electrode and oxygen (or another oxygen-containing product such as hydrogen peroxide) will be produced at the counter electrode, so it may be desirable to add nitrogen gas at the working electrode and remove oxygen gas at the counter electrode; alternative the nitrogen may simply be suppled from ambient air, and oxygen released to the ambient air. Water may be supplied as a gas by humidifying the supplied nitrogen, or by addition of liquid water to the electrolyte solution.

The iron-containing Chevrel phase material is a very active and efficient catalyst for the electrochemical production of ammonia. The mass-normalized yield rate of NH₃ production in the electrochemical production of ammonia from water and nitrogen gas is preferably at least 25 μg/h per mg of catalyst, more preferably 50 μg/h per mg of catalyst, and most preferably at least 100 μg/h per mg of catalyst including 70-100 μg/h per mg of catalyst. The Faradaic efficiency of the reaction in the electrochemical production of ammonia from water and nitrogen gas is preferably at least 12%, more preferably at least 20% and most preferably at least 25%, including 12% to 25%. Greater values of mass-normalized yield rate and Faradaic efficiency may be achieved by adjust the temperature and/or pressure of the reaction, as well as amount of nitrogen dissolved in the electrolyte solution.

As illustrated in FIG. 4, a facility for manufacturing ammonia, 50, includes a system for producing ammonia electrochemically, 52, a water feed, 54, an optional nitrogen feed, 56, an ammonia dispenser for dispensing ammonia, ammonia containing gas, or ammonia or ammonium containing solution, 58, and an optional oxygen dispenser, 60. Also illustrated are power sources from providing electrical power to the facility, which may include one or more sources including hydroelectric or geothermal power source, 62, solar power source, 64, wind power source, 66, and municipal or generator power source, 68, which are electrically connected to the facility by electrical cables, 70.

EXAMPLES Example 1

Materials and Methods

Chemicals: Iron sulfide (99%, Sigma-Aldrich), molybdenum disulfide (>99%, Sigma-Aldrich), molybdenum (≥95%, VWR), copper sulfide (>99%, Alfa Aesar), nitric acid (98%, Fisher Scientific), hydrochloric acid (HCl, 37.3% Fisher Scientific), ammonium hydroxide (30%, Sigma-Aldrich), ammonium chloride (≥98%, Sigma-Aldrich), sodium hypochlorite solution (NaClO, available chlorine>5.0%), sodium hydroxide (98.8%, Fisher Scientific), sodium citrate dehydrate (>99%, Fisher Scientific), phenol (>99%, ACROS), sodium nitroferricyanide dihydrate (>99%, ACROS), p-dimethylaminobenzaldehyde (>99.9%, Fisher Scientific), and NAFION™ (5 wt %, Sigma-Aldrich). All chemicals were used as-received without further purification.

Synthesis Fe₂Mo₆S₈. In a typical procedure, 0.88 g FeS, 1.44 g Mo and 2.40 g MoS₂ were loaded to a ball milling jar inside a Ar-filled glove box. The jar was then sealed and transferred out and was milled with a high-energy mechanical miller (8000 M, SPEX SamplePrep, USA) for 9 hours. After this process, the powders were collected and transferred to a Lindberg Blue tube furnace. The tube furnace was heated to 1000° C. for 10 hours under the flow of 100 sccm Ar, and was then cooled to room temperature. The product was Fe₂Mo₆S₈.

Preparation of Mo₆S₈: Mo₆S₈ particles were prepared from Cu₂Mo₆S₈ by acid leaching off Cu²⁺. The Cu₂Mo₆S₈ powders were synthesized using the method as detailed in our previous publications.⁴⁷ In a typical procedure, 3.0 g of the as-synthesized Cu₂Mo₆S₈ powders were added to 50 ml of the mixed acid containing 8.0 M HCl and 0.05 M HNO₃. The mixture was stirred at room temperature for 1 day and was then collected by centrifuge, washed repeatedly with water until pH neutral, and was then dried under vacuum at 100° C.

Materials Characterization

Powder XRD was recorded using a Rigaku Miniflex diffractometer, with a Cu Kα radiation (λ=1.5406 Å, 30 kV, 15 mA). The microstructures of the prepared samples were collected by SEM images on a field-emission Hitachi S-4700-II microscope. TEM images were obtained on a JEOL JEM-2100F microscope at 200 kV. XPS measurements were performed on a Thermal Escalab 250 X-ray photoelectron spectrometer, and the binding energies were calibrated by assigning the C 1s peak at 284.5 eV. Raman spectra were collected on a Renishaw inVia Raman microscope with a 532 nm laser. Specific surface area and pore size distribution in powder samples were analyzed using a Micromeritics Tristar 3000 analyzer. UV-vis absorbance spectra were measured applying a SHIMADZU UV-1800 UV-vis spectrophotometer.⁵⁷Fe Mössbauer spectra were recorded in transmission mode using a ⁵⁷Co/Rh radiation source with a Kr proportional counter on a constant acceleration drive (SEECO, Edina, Minn.). The isomer shifts were given relative to α-Fe and the velocity was calibrated using a α-Fe foil. The Fe₂Mo₆S₈ sample was prepared by pressing ca. 200 mg powders into pellets (approximately 15 mm in diameter and 2 mm in thickness), and were measured at ambient temperature with a velocity range of ±8.0 mm s⁻¹.

Electrochemical Studies

All electrochemical tests were performed using a Pine bipotentiostat electrochemical workstation (Pine Instruments, USA) and a two-compartment H-cell separated by a NAFION™ 211 membrane (Fuelcell Store) at room temperature. A standard three-electrode configuration was used, with a saturated calomel reference electrode (SCE) and a rotating disk glassy carbon electrode coated with relevant catalysts working electrode in one compartment of the H-cell and the a graphite rod counter electrode was in the other compartment of the H-cell. The electrolytes used in this work were 0.5 M Na₂SO₄ and 0.1 M sodium citrate buffer (pH=4.0). The catalyst ink was prepared by mixing 10 mg catalyst with 1.9 mL water/isopropanol (1:3 v/v) and 0.1 ml 5 wt % NAFION™ solution. The mixture was sonicated for 30 min to form uniform dispersion. To prepare the working electrode, 8 μL of the catalyst ink was drop-casted onto a 5 mm glassy carbon rotating disk electrode. The mass loading was 0.2 mg cm⁻². The linear sweep voltammograms (LSVs) were collected at 5 mV s⁻¹. All LSV curves were iR corrected using the standard equation of E_(corrected)=E_(measured)−iR_(s), where E_(corrected), E_(measured) and i refers to the iR-corrected potential, measured potential and recorded current, respectively. R_(s) is the equivalent resistance measured by electrochemical impedance spectroscopy (EIS) employing the same electrode configuration. The potentiostat tests were performed at different potentials including −0.15, −0.20, −0.25, −0.30, −0.40, −0.50 vs. RHE. The electrochemically active surface area, A_(ECSA), was estimated using the equation of A_(ECSA)=C_(d)/C_(s), where the C_(dl) and C_(s) refer to the double layer capacitance and a specific capacitance value of 22 μF cm⁻², respectively.^(48, 49)

Ammonia detection: The indophenol blue method was employed to quantify the concentration of ammonia generated in the electrochemical cell.⁵⁰ Three solutions were prepared for this method, including the phenol solution, prepared by adding 20.0 g phenol in 100 ml ethanol; the sodium nitroprusside solution, prepared by dissolving 1.0 g sodium nitroferricyanide in 200 ml water in an amber bottle; the alkaline reagent, prepared by dissolving 100 g sodium citrate and 5 g sodium hydroxide in 500 ml water. The oxidizing solution was prepared fresh each time by mixing 10 ml stock alkaline reagent and 2.5 ml sodium hyphchlorite (NaClO) solution. The calibration curve was established using standard NH₄Cl stock solution, in the concentration range of 0, 0.5, 1, 5, 10, 50 and 100 μM. In each test 3.0 ml of the standard solution with the designated NH₄Cl concentration was mixed with 0.12 ml phenol solution, 0.12 ml sodium nitroprusside solution and 0.3 ml oxidizing solution. The mixed solution was mixed using a vortex mixer and allowed to stand at room temperature for at least 1 hour in the dark. The concentration of indophenol blue was determined using the absorbance at the wavelength of 640 nm. The fitting curve shows good linear relationships. The ammonia concentration in the electrolyte was estimated using the calibration curve, and samples were prepared by mixing 3.0 ml electrolyte, 0.12 ml phenol solution, 0.12 ml sodium nitroprusside solution and 0.3 ml oxidizing solution. The same aging time and analysis protocol were used.

¹⁵N isotopic labeling experiment. ¹⁵N₂ was used as the feeding gas, a low velocity gas. After electrolysis at −0.25 V vs. RHE for 10 hours. 0.9 ml of the concentrated electrolyte was collected, followed by adding 0.1 ml of D₂O as the internal standard. The produced ammonia was quantified using ¹H nuclear magnetic resonance measurements (¹H NMR; Bruker 300 MHz).

N₂H₄ detection: The concentration of the hydrazine was determined by the Watt and Chrisp method.⁵¹ The coloring reagent was prepared by mixing 300 ml ethanol, 30 ml concentrated HCl, and 5.99 g p-dimethylaminobenzaldehyde. In a typical test, 5 ml of the electrolyte solution after the NRR potentiostatic test was collected and mixed with 5 ml of the coloring reagent solution. The resulting solution was stirred for 10 min, and its absorbance was measured at a wavelength of 455 nm. Standard hydrazine monohydrate solutions at a series of concentrations in 0.5 M Na₂SO₄ and 0.1 M sodium citrate buffer were prepared to establish the calibration curve.

Faraday efficiency: The Faraday efficiency (FE) and mass-normalized yield rate of NH₃ production were calculated using the following equation:⁵²

FE(NH₃)=[3F×c(NH₃)×V]/Q

Yield Rate_(mass)(NH₃)=[17×c(NH₃)×V]/(t×m)

where F is the Faraday constant (96485 C mol⁻¹), t is the electrolysis time, m is the loading mass of the catalyst, Q is the total charge passed through the electrode, V is the volume of the electrolyte and c(NH₃) is the quantified ammonia concentration.

The surface-area normalized yield rate of NH₃ production was quantified using the equation:

Yield Rate_(ESCA)(NH₃)=[17×c(NH₃)×V]/(t×A_(ESCA))

where A_(ESCA) is the electrochemically active surface area.

Results and Discussion

Phase pure Fe₂Mo₆S₈ catalysts were synthesized using a two-step method. First, FeS, Mo and MoS₂ powders with the molar ratio of 2:3:3 were loaded in a ball-milling container inside an Ar-filled glove box and were milled for 9 hours using a high-energy mechanical miller. The milled mixture was collected and transferred to a tube furnace and calcined under Ar at 1000° C. for 10 hours, which resulted in formation of nearly pure Fe₂Mo₆S₈ powders (typically ˜4.0 g). The powder X-ray diffraction analysis of the as-synthesized particles revealed patterns that index well with the rhombohedral structure of the Chevrel phase with the R3 space group (FIG. 5).³² The structure can be viewed as three-dimensional frameworks with interconnected quasi-rigid Mo₆S₈ clusters, each of which has a Mob octahedral surrounded by a Sa cube. The Fe-ions distribute among sites A and B in channels structured by Mo₆S₆ clusters as illustrated in FIG. 1, which is comparable to Cu-ions in Cu₂Mo₆S₈ from previous studies given their similar cation size.³³ The scanning (SEM) and transmission electron, microscope (TEM) images revealed that the particles have a wide size distribution of 10˜100 nm (FIG. 6, FIG. 7 and FIG. 8). Most particles were aggregated due to the nature of ball-milling based synthesis and have rough surfaces with abundant edges and corners, which usually have higher electrocatalytic activity compared with flat basal planes.³⁴ The high resolution TEM images revealed a large (101) rhombohedral lattice spacing of 6.67 Å, which is slightly larger than the 6.45 Å of Mo₆S₈ due to Fe-ions induced lattice expansions (FIG. 8).³⁵ in addition, the produced Fe₂Mo₆S₈ catalysts exhibited a Brunauer-Emmett-Teller (BET) specific surface area of only ˜2 m² g⁻¹ due to their relatively large size compared with conventional catalysts (FIG. 9).

The produced Fe₂Mo₆S₈ particles were directly studied without further treatments as electrocatalysts for ambient N₂ reduction reaction (NRR). A two compartment H-type electrochemical cell separated by a piece of NAFION™ 211 membrane was employed together with a three-electrode setup to avoid contaminations from electrolysis products formed at the counter electrode in the analysis of N₂ conversion products.¹⁵ A saturated calomel reference electrode (SCE) and a 4.0 mm rotating disk glassy carbon working electrode (RDE) were mounted in one compartment and a graphite rod counter electrode was mounted in the other compartment of the H-cell. The electrolyte was an optimized aqueous solution of 0.5 M Na₂SO₄ mixed with 0.1 M sodium citrate (pH˜4.0), where the Na₂SO₄ was added to ensure ionic conductivity. The catalyst loading was 0.2 mg cm⁻² for all experiments and no carbon additive was employed in the formulation of the catalyst inks in order to assess the true activities of Fe₂Mo₆S₈ and avoid interferences from carbon.³⁶ In addition, all electrolysis was tested with 800 rpm to ensure uniform N₂ feed to electrocatalysts.

The catalysts were first activated in the Ar-saturated electrolyte using cyclic voltammetry between −0.1 and 0.2 V (vs. RHE, at 10 mV s⁻¹) until the voltammogram stabilized. This small voltage range is essential to avoid anodic leaching of Fe-ions from Fe₂Mo₆S₈. FIG. 10 compares the iR-corrected anodeic linear sweep voltammograms (LSV) in Ar and N₂ saturated electrolytes at 5 mV s⁻¹. The voltammogram under N₂ exhibited noticeably larger current densities, suggesting Faradaic contributions from electrochemical reactions that are associated with N₂. FIG. 11 compares Tafel plots derived from the LSV results. The voltammogram under Ar exhibited a Tafel slope of 203 mV dec⁻¹ for the HER, which suggest sluggish HER kinetics that is beneficial for improving NH₃ selectivity.²⁷ In contrast, the slope decreased to 170 mV dec⁻¹ in the N₂-saturated electrolyte, clearly suggesting electrochemical processes that N₂-assisted H_(ad)* generation and/or removal from the catalyst surface.

We then performed chronoamperometry tests under N₂ bubbling and analyzed the collected electrolytes to validate NH₃ synthesis activity of the Fe₂Mo₆S₈ catalyst. FIG. 12 presents the chronoamperometric curves for 2 hours under the applied potentials of −0.15, −0.20, −0.25, −0.30, −0.40 (vs. RHE). Notably, the current density at different potentials remained relatively stable throughout all testing, revealing good stability of Fe₂Mo₆S₈ and can be ascribed to its unique framework structure.³⁷ The collected electrolytes were analyzed for NH₃ concentration using the indophenol blue method in reference to standard NH₄Cl solutions (FIG. 13 and FIG. 14),³⁸ and the results are presented in FIG. 15 for NH₃ yield rates and Faradaic efficiencies (FE) under different potentials (with the curve being Faradaic efficiency, and the bars showing the yield in the figure). The Fe₂Mo₆S₈ catalyst exhibited remarkable activities for catalyzing N₂ to NH₃ despite its orders-of-magnitude lower specific surface area compared with prevalent electrocatalysts.³⁹ The formation of ammonia was clearly detected at −0.15 V (vs. RHE), with an FE of 10.3% and NH₃ production rate of 17 μg h⁻¹ mg_(cat) ⁻¹. The electrolysis at −0.20 V vs. RHE, on the other hand, exhibited the best combination of FE and NH₃ production rate. While the FE increased slightly to 12.5%, the NH₃ yield rate increased approximately five-fold to ˜70 μg h⁻¹ mg_(cat) ⁻¹ and suggesting boosted N₂ conversion despite the current density only increasing two-fold to ˜0.25 mA cm⁻². This small current density is vital for high FE and was obtained by inhibiting the HER through the synergistic Fe—Mo₆S₈ modulations as discussed below. As expected, the FE decreased gradually at more negative potentials, reaching to ˜5% at −0.25 V and ˜3% at −0.30 V due to rapid rising of HER activity. Interestingly, the NH₃ production rate remained relatively stable at 50˜70 μg h⁻¹ mg_(cat) ⁻¹. Similarly stable production rates were also observed with the Pd/C catalyst suggesting the production rate was likely limited by soluble N₂ in the electrolyte.¹⁵ In addition, we performed control experiments with Ar-saturated electrolyte or N₂-saturated electrolyte but without Fe₂Mo₆S₈, and confirmed that no detectable NH₃ was observed in both cases.

FIG. 16 presents the relationship between ammonia concentration and electrolysis duration at the optimal potential of −0.20 V vs. RHE. A nearly linear relationship was observed, with a slope corresponding to the yielding rate of ˜70 μg h⁻¹ mg_(cat) ⁻¹ at this potential (dashed line), further demonstrating that Fe₂Mo₆S₈ is a highly stable catalyst for ambient electrochemical N₂ conversion. A slight departure from the linear relationship was observed after 20 hours, this is likely due to the gradual pH change of the electrolyte solution (from 4.0 to 4.9 after 50 hours of electrolysis), and this changed the activity of H₃O⁺ causing the slight departure from the ideal condition. We further employed the ¹⁵N₂ isotope labeling experiment to confirm the observed NH₃ was indeed generated from N₂ (FIG. 17).⁴⁰ The ¹H NMR spectrum of electrolyte collected with ¹⁵N₂ feeding gas revealed a doublet coupling that corresponds to ¹⁵NH₄ ⁺. In comparison, the triplet coupling for¹⁴NH₄ ⁺ was observed in the electrolyte produced with ¹⁴N₂ feeding gas. This comparison clearly confirms the true activity of Fe₂Mo₆S₈ catalysts for N₂ to NH₃ conversion. In addition, X-ray photoelectron spectroscopy (XPS) confirmed the Fe₂Mo₆S₈ catalysts did not contain any nitrogen-species, and therefore the formation of ammonia from catalyst decomposition is unlikely to be a concern (FIG. 18).⁴¹ Besides NH₃, the N₂ conversion could possibly proceed along an alternative pathway and produce N₂H₄. In this study, no detectable formation of N₂H₄ was recorded using the Watt and Chrisp method (FIG. 19).⁴² The UV-Vis absorption spectra of samples prepared from fresh electrolyte and electrolyte after 2 hours of electrolysis at −0.2 V vs RHE did not exhibit any features corresponding to N₂H₄, confirming nearly exclusive electrocatalytic N₂ conversion to NH₃.

The Fe₂Mo₆S₈ catalyst is a pure compound and its active sites should be based on its multi-element Fe, Mo and S configuration, which enables superior N₂ conversion (FIG. 20). Indeed, previous theoretical studies outlined that Fe-ions moderate the catalytic activities of Mo₆S₈ via three effects: the ligand effect where Fe donates electrons to Mo₆S₈ and increases electron densities on Mo and S, the ensemble effect where Fe directly participates in binding key reactants/intermediates, and the confinement effect where the unique configuration of Fe, Mo and S spatially confine reactions to proceed along specific pathways.^(37, 43) Here, the Fe—Mo₆S₈ modulations were studied by spectroscopic characterizations of Fe₂Mo₆S₈, and the spectra from Mo₆S₈ were collected as references (FIG. 21). FIG. 22 presents the Raman spectrum of Fe₂Mo₆S₈, which when compared to the spectra of FeS and Mo₆S₈ confirms resonance peaks from the E_(g) and A_(g) vibration modes of Fe, S and Mo.⁴⁴ The type of Fe species present in Fe₂Mo₆S₈ were established using Mössbauer characterizations at room temperature (FIG. 23). The fitting of the spectrum suggests three major types of Fe species (see Table 1 for detailed fitting parameters), with 21.3% as FeMo₁ and FeMo₀ alloys (0 and 1 refer to number of nearby Mo atom), 33.6% as Fe₀ and 45.1% as Fe ions that can be described as Fe²⁺/Fe³⁺ mixtures.^(45, 46) The presence of multiple Fe species can attribute to the unique versatile framework structure of the Chevrel phase structure. The exact roles of each Fe species on N₂ conversion are unclear and will be studied in following works. We do note, however, that similar Fe species were described in the FeMo cofactor protein and Fe-ions with multiple oxidation states are probably essential for synergistic N₂ absorption and activation.^(9, 30) The high-resolution XPS spectra revealed increased electron densities on Mo and S in Fe₂Mo₆S₈ and therefore confirmed the ligand role of Fe on modulating Mo₆S₈ (FIG. 24 and FIG. 25). The results suggest S in Fe₂Mo₆S₈ has 0.5˜0.8 eV lower binding energies than in Mo₆S₈; for Mo, a portion of Mo³⁺ in Mo₆S₈ was reduced to Mo²⁺ in Fe₂Mo₆S₈. These results agree well with theoretical predictions, and the ligand effect play keys roles for inhibiting HER and improving NH₃ selectivity.

TABLE 1 Fitting parameters for Fe-species identified by Mossbauer spectrum (FIG. 23)^(53,54) Isomer Quadrupole Magnetic field, Component shift, mm/s splitting, mm/s KOe Area, % FeMo₁ 0.023 0.130 295.0 8.43 FeMo₀ 0.001 0.043 332.4 12.87 Fe²⁺/Fe³⁺ 0.746 1.703 0 45.09 Fe⁰ 0.014 0.001 0 33.60

The ensemble effect of Fe-species in Fe₂Mo₆S₈ plays a vital role in N₂ to NH₃ conversion,⁴³ as control catalysts with the same Chevrel phase structure but different composition, including Cu₂Mo₆S₈, Mn₂Mo₆S₈ (FIG. 26 and FIG. 27 for XRD) and Mo₆S₈, all exhibited much lower activity at −0.2V vs. RHE (FIG. 28). This comparison suggests that Fe must participate in the absorption and/or conversion of N₂ while Mo/S atoms alone were unable to activate N₂. Rather, they assist generation and stabilization of the hydrogen intermediate H_(ad)* by providing strong S—H binding, which is essential to inhibit the undesired HER.³⁴ This inhibiting effect was strengthened with metal cations in Mo₆S₈ (the ligand effect), as higher overpotentials and smaller current densities for the HER were observed for Fe₂Mo₆S₈, Cu₂Mo₆S₈ and Mn₂Mo₆S₈ (FIG. 29). Notably, Fe-ions exhibited the most pronounced ligand effect with the Fe₂Mo₆S₈ catalysts having the most sluggish kinetics for the HER, with a ˜100 mV higher overpotential and ˜50% lower current density compared with Mo₆S₈, along with a much larger Tafel slope (203 vs 114 mV dec⁻¹ for Mo₆S₈, FIG. 11 and FIG. 30). The ligand effect of Fe increases electron density on S as discussed above, which strengthens the S—H_(ad)* binding and extends the lifespan of H_(ad)* for robust association with *N₂. Overall, these understandings could be rationalized into a tentative mechanistic view of the reaction pathway. In this model, the highly conductive Fe₂Mo₆S₈ functions as a true multi-site catalyst and has dedicated binding sites for *N₂ and *H on Fe/Mo and S, respectively. The unique spatial geometry of Fe, Mo and S in Fe₂Mo₆S₈ catalysts confine these two key intermediates in close proximity for promoted association and formation of NH₃, leading to a remarkable NH₃ production rate of 70 μg h⁻¹ mg_(cat) ⁻¹. Furthermore, this rate translates into a surface normalized production rate of 3.5 μg h⁻¹ cm_(cat) ⁻² as the catalyst has a low BET surface area of ˜2.0 m² g⁻¹, which is orders of magnitude higher than the typical existing catalysts. The production rate could be further enhanced by advanced materials synthesis for nanoscale Fe₂Mo₆S₈, and therefore could provide some new insights on the design principles of selective catalysts for N₂ to NH₃ conversion.

Conclusion

In summary, we describe outstanding activities of the Iron-containing Chevrel phase chalcogenides, such as Fe₂Mo₆S₈, for highly selective electrochemical conversion of N₂ to NH₃ under ambient condition. The activities arise from the unique structure of these materials, which provides true multi-active binding sites for separate binding and activating key precursor molecules, including the Fe/Mo-sites for activating N₂ and the S-site for binding with H and inhibiting the undesired HER. The geometry of Fe, Mo and S spatially confines these intermediates to ensure facile hydrogenation of *N₂ for promoted NH₃ formation, reaching a high NH₃ production rate of 70 μg h⁻¹ mg_(cat) ⁻¹ that translate into a remarkable surface area normalized rate of 3.5 μg h⁻¹ cm_(cat) ⁻². The observation of outstanding activity with low surface area Fe₂Mo₆S₈ catalysts is surprising, and suggests the catalyst is intrinsically active for N₂ conversion.

Example 2: Nitrogen Reduction to Ammonia on Fe₄Mo₆S₈

Synthesis:

Phase pure Fe₄Mo₆S₈ catalysts were synthesized using a two-step method. First, FeS, Mo and MoS₂ powders with the molar ratio of 2:2:1 (typically 1.76 g FeS, 1.92 g Mo and 1.5 g MoS₂) were loaded in a ball-milling container inside an Ar-filled glove box and were milled for 9 hours using a high-energy mechanical miller. The milled mixture was collected and transferred to a tube furnace and calcined under Ar at 1000° C. for 10 hours, which resulted in formation of nearly pure Fe₄Mo₆S₈ powders. The X-ray powder diffraction pattern of Fe₄Mo₆S₈ is shown in FIG. 32, and SEM images are shown in FIG. 33 and FIG. 34. Fe_(x)Mo₆S₈, where x=1 and 6 were prepared in a similar fashion, but some impurities were present in the product.

Fe₄Mo₆S₈ was tested as a catalyst for the nitrogen reduction of ammonia, using the same testing protocols as described in Example 1. Linear sweep voltammetry (LSV) of Fe₄Mo₆S₈ in Ar and N₂ saturated electrolytes are shown in FIG. 35. The larger cathodic reduction current suggests Faradaic reactions associated with nitrogen. Comparison of Faradaic efficiencies for ammonia production and average ammonia yield at different applied currents is shown in FIG. 36. The highest Faradaic efficiency was 25% and the highest ammonia yield was ˜102 μg h⁻¹ mg_(cat) ⁻¹. The Fe₄Mo₆S₈ catalyst has good stability for continuous production of ammonia. No appreciable decay in Faradaic efficiency and average ammonia yield was observed within 10 cycles, as shown in FIG. 37.

REFERENCES

-   1. Erisman, J. W.; Sutton, M. A.; Galloway, J.; Klimont, Z.;     Winiwarter, W., How a century of ammonia synthesis changed the     world. Nat Geosci 2008, 1 (10), 636-639. -   2. Guo, J.; Chen, P., Catalyst: NH₃ as an Energy Carrier. Chem 2017,     3 (5), 709-712. -   3. Wang, L.; Xia, M.; Wang, H.; Huang, K.; Qian, C.; Maravelias, C.     T.; Ozin, G. A., Greening Ammonia toward the Solar Ammonia Refinery.     Joule 2018, 2 (6), 1055-1074. -   4. Kandemir, T.; Schuster, M. E.; Senyshyn, A.; Behrens, M.;     Schlögl, R., The Haber-Bosch Process Revisited: On the Real     Structure and Stability of “Ammonia Iron” under Working Conditions.     Angewandte Chemie International Edition 2013, 52 (48), 12723-12726. -   5. Soloveichik, G., Electrochemical synthesis of ammonia as a     potential alternative to the Haber-Bosch process. Nat Catal 2019, 2     (5), 377-380. -   6. Smith, C.; Hill, A. K.; Torrente-Murciano, L., Current and future     role of Haber-Bosch ammonia in a carbon-free energy landscape.     Energy & Environmental Science 2020, 13 (2), 331-344. -   7. Zheng, J.; Lyu, Y.; Qiao, M.; Wang, R.; Zhou, Y.; Li, H.; Chen,     C.; Li, Y.; Zhou, H.; Jiang, S. P.; Wang, S., Photoelectrochemical     Synthesis of Ammonia on the Aerophilic-Hydrophilic Heterostructure     with 37.8% Efficiency. Chem 2019, 5 (3), 617-633. -   8. Lazouski, N.; Schiffer, Z. J.; Williams, K.; Manthiram, K.,     Understanding Continuous Lithium-Mediated Electrochemical Nitrogen     Reduction. Joule 2019, 3 (4), 1127-1139. -   9. Foster, S. L.; Bakovic, S. I. P.; Duda, R. D.; Maheshwari, S.;     Milton, R. D.; Minteer, S. D.; Janik, M. J.; Renner, J. N.;     Greenlee, L. F., Catalysts for nitrogen reduction to ammonia. Nat     Catal 2018, 1 (7), 490-500. -   10. Rosca, V.; Duca, M.; de Groot, M. T.; Koper, M. T. M., Nitrogen     Cycle Electrocatalysis. Chemical Reviews 2009, 109 (6), 2209-2244. -   11. Hu, L.; Xing, Z.; Feng, X. F., Understanding the     Electrocatalytic Interface for Ambient Ammonia Synthesis. Acs Energy     Letters 2020, 5 (2), 430-436. -   12. Singh, A. R.; Rohr, B. A.; Statt, M. J.; Schwalbe, J. A.;     Cargnello, M.; Nørskov, J. K., Strategies toward Selective     Electrochemical Ammonia Synthesis. ACS Catalysis 2019, 9 (9),     8316-8324. -   13. Hawtof, R.; Ghosh, S.; Guarr, E.; Xu, C.; Mohan Sankaran, R.;     Renner, J. N., Catalyst-free, highly selective synthesis of ammonia     from nitrogen and water by a plasma electrolytic system. Sci Adv     2019, 5 (1), eaat5778. -   14. Hu, L.; Khaniya, A.; Wang, J.; Chen, G.; Kaden, W. E.; Feng, X.,     Ambient Electrochemical Ammonia Synthesis with High Selectivity on     Fe/Fe Oxide Catalyst. ACS Catalysis 2018, 8 (10), 9312-9319. -   15. Wang, J.; Yu, L.; Hu, L.; Chen, G.; Xin, H.; Feng, X., Ambient     ammonia synthesis via palladium-catalyzed electrohydrogenation of     dinitrogen at low overpotential. Nat Commun 2018, 9 (1), 1795. -   16. Li, L.; Tang, C.; Xia, B.; Jin, H.; Zheng, Y.; Qiao, S.-Z.,     Two-Dimensional Mosaic Bismuth Nanosheets for Highly Selective     Ambient Electrocatalytic Nitrogen Reduction. ACS Catalysis 2019, 9     (4), 2902-2908. -   17. Patil, S. B.; Chou, H.-L.; Chen, Y.-M.; Hsieh, S.-H.; Chen,     C.-H.; Chang, C.-C.; Li, S.-R.; Lee, Y.-C.; Lin, Y.-S.; Li, H.;     Chang, Y. J.; Lai, Y.-H.; Wang, D.-Y., Enhanced N₂ affinity of     1T-MoS₂ with a unique pseudo-six-membered ring consisting of     N—Li—S—Mo—S—Mo for high ambient ammonia electrosynthesis     performance. Journal of Materials Chemistry A 2021. -   18. Tao, H.; Choi, C.; Ding, L.-X.; Jiang, Z.; Han, Z.; Jia, M.;     Fan, Q.; Gao, Y.; Wang, H.; Robertson, A. W.; Hong, S.; Jung, Y.;     Liu, S.; Sun, Z., Nitrogen Fixation by Ru Single-Atom     Electrocatalytic Reduction. Chem 2019, 5 (1), 204-214. -   19. Pang, F.; Wang, Z.; Zhang, K.; He, J.; Zhang, W.; Guo, C.; Ding,     Y., Bimodal nanoporous Pd₃Cu₁ alloy with restrained hydrogen     evolution for stable and high yield electrochemical nitrogen     reduction. Nano Energy 2019, 58, 834-841. -   20. Han, L.; Liu, X.; Chen, J.; Lin, R.; Liu, H.; Lu, F.; Bak, S.;     Liang, Z.; Zhao, S.; Stavitski, E.; Luo, J.; Adzic, R. R.; Xin, H.     L., Atomically Dispersed Molybdenum Catalysts for Efficient Ambient     Nitrogen Fixation. Angew Chem Int Ed Engl 2019, 58 (8), 2321-2325. -   21. Deng, J.; Iñiguez, J. A.; Liu, C., Electrocatalytic Nitrogen     Reduction at Low Temperature. Joule 2018, 2 (5), 846-856. -   22. Zhang, L.; Ji, X.; Ren, X.; Ma, Y.; Shi, X.; Tian, Z.; Asiri, A.     M.; Chen, L.; Tang, B.; Sun, X., Electrochemical Ammonia Synthesis     via Nitrogen Reduction Reaction on a MoS₂ Catalyst: Theoretical and     Experimental Studies. Adv Mater 2018, 30 (28), e1800191. -   23. Cheng, H.; Ding, L. X.; Chen, G. F.; Zhang, L.; Xue, J.; Wang,     H., Molybdenum Carbide Nanodots Enable Efficient Electrocatalytic     Nitrogen Fixation under Ambient Conditions. Adv Mater 2018, 30 (46),     e1803694. -   24. Chen, G. F.; Cao, X.; Wu, S.; Zeng, X.; Ding, L. X.; Zhu, M.;     Wang, H., Ammonia Electrosynthesis with High Selectivity under     Ambient Conditions via a Li(+) Incorporation Strategy. J Am Chem Soc     2017, 139 (29), 9771-9774. -   25. Bao, D.; Zhang, Q.; Meng, F. L.; Zhong, H. X.; Shi, M. M.;     Zhang, Y.; Yan, J. M.; Jiang, Q.; Zhang, X. B., Electrochemical     Reduction of N₂ under Ambient Conditions for Artificial N₂ Fixation     and Renewable Energy Storage Using N₂/NH₃ Cycle. Adv Mater 2017, 29     (3). -   26. Singh, A. R.; Rohr, B. A.; Schwalbe, J. A.; Cargnello, M.; Chan,     K.; Jaramillo, T. F.; Chorkendorff, I.; Nørskov, J. K.,     Electrochemical Ammonia Synthesis—The Selectivity Challenge. ACS     Catalysis 2017, 7 (1), 706-709. -   27. Lu, K.; Liu, Y.; Lin, F.; Cordova, I. A.; Gao, S.; Li, B.; Peng,     B.; Xu, H.; Kaelin, J.; Coliz, D.; Wang, C.; Shao, Y.; Cheng, Y.,     LixNiO/Ni Heterostructure with Strong Basic Lattice Oxygen Enables     Electrocatalytic Hydrogen Evolution with Pt-like Activity. Journal     of the American Chemical Society 2020, 142 (29), 12613-12619. -   28. Suryanto, B. H. R.; Wang, D.; Azofra, L. M.; Harb, M.; Cavallo,     L.; Jalili, R.; Mitchell, D. R. G.; Chatti, M.; MacFarlane, D. R.,     MoS₂ Polymorphic Engineering Enhances Selectivity in the     Electrochemical Reduction of Nitrogen to Ammonia. ACS Energy Letters     2018, 4 (2), 430-435. -   29. Burgess, B. K.; Lowe, D. J., Mechanism of Molybdenum     Nitrogenase. Chem Rev 1996, 96 (7), 2983-3012. -   30. Buscagan, T. M.; Rees, D. C., Rethinking the Nitrogenase     Mechanism: Activating the Active Site. Joule 2019, 3 (11),     2662-2678. -   31. Banerjee, A.; Yuhas, B. D.; Margulies, E. A.; Zhang, Y.; Shim,     Y.; Wasielewski, M. R.; Kanatzidis, M. G., Photochemical nitrogen     conversion to ammonia in ambient conditions with FeMoS-chalcogels. J     Am Chem Soc 2015, 137 (5), 2030-4. -   32. Mao, M.; Lin, Z.; Tong, Y.; Yue, J.; Zhao, C.; Lu, J.; Zhang,     Q.; Gu, L.; Suo, L.; Hu, Y.-S.; Li, H.; Huang, X.; Chen, L., Iodine     Vapor Transport-Triggered Preferential Growth of Chevrel Mo₆S₈     Nanosheets for Advanced Multivalent Batteries. ACS Nano 2020, 14     (1), 1102-1110. -   33. Levi, E.; Gershinsky, G.; Aurbach, D.; Isnard, O.,     Crystallography of Chevrel phases, MMo₆T₈ (M=Cd, Na, Mn, and Zn,     T=S, Se) and their cation mobility. Inorg Chem 2009, 48 (18),     8751-8. -   34. Liu, M.; Hybertsen, M. S.; Wu, Q., A Physical Model for     Understanding the Activation of MoS₂ Basal-Plane Sulfur Atoms for     the Hydrogen Evolution Reaction. Angewandte Chemie International     Edition 2020, 59 (35), 14835-14841. -   35. Lu, K.; Gao, S.; Li, G.; Kaelin, J.; Zhang, Z.; Cheng, Y.,     Regulating Interfacial Na-Ion Flux via Artificial Layers with Fast     Ionic Conductivity for Stable and High-Rate Na Metal Batteries. ACS     Materials Letters 2019, 1 (3), 303-309. -   36. Sheng, H.; Hermes, E. D.; Yang, X.; Ying, D.; Janes, A. N.; Li,     W.; Schmidt, J. R.; Jin, S., Electrocatalytic Production of H₂O₂ by     Selective Oxygen Reduction Using Earth-Abundant Cobalt Pyrite     (CoS₂). ACS Catalysis 2019, 9 (9), 8433-8442. -   37. Zhang, H. T.; Liu, C.; Liu, P.; Hu, Y. H., Mo₆S₈-based     single-metal-atom catalysts for direct methane to methanol     conversion. J Chem Phys 2019, 151 (2), 024304. -   38. Bolleter, W. T.; Bushman, C. J.; Tidwell, P. W.,     Spectrophotometric Determination of Ammonia as Indophenol. Anal.     Chem. 1961, 33, 592-594. -   39. Zhao, S.; Lu, X.; Wang, L.; Gale, J.; Amal, R., Carbon-Based     Metal-Free Catalysts for Electrocatalytic Reduction of Nitrogen for     Synthesis of Ammonia at Ambient Conditions. Adv Mater 2019, 31(13),     e1805367. -   40. Andersen, S. Z.; Colic, V.; Yang, S.; Schwalbe, J. A.;     Nielander, A. C.; McEnaney, J. M.; Enemark-Rasmussen, K.; Baker, J.     G.; Singh, A. R.; Rohr, B. A.; Staff, M. J.; Blair, S. J.;     Mezzavilla, S.; Kibsgaard, J.; Vesborg, P. C. K.; Cargnello, M.;     Bent, S. F.; Jaramillo, T. F.; Stephens, I. E. L.; Norskov, J. K.;     Chorkendorff, I., A rigorous electrochemical ammonia synthesis     protocol with quantitative isotope measurements. Nature 2019, 570     (7762), 504-508. -   41. Hu, B.; Hu, M.; Seefeldt, L.; Liu, T. L., Electrochemical     Dinitrogen Reduction to Ammonia by Mo₂N: Catalysis or Decomposition?     ACS Energy Letters 2019, 4 (5), 1053-1054. -   42. Watt, G. W.; Chrisp, J. D., Spectrophotometric Method for     Determination of Hydrazine. Anal Chem 1952, 24 (12), 2006-2008. -   43. Liu, C.; Liu, P., Mechanistic Study of Methanol Synthesis from     CO₂ and H₂ on a Modified Model Mo₆S₈ Cluster. ACS Catalysis 2015, 5     (2), 1004-1012. -   44. Holmgren, D. J.; Demers, R. T.; Klein, M. V.; Ginsberg, D. M.,     Raman study of phonons in Chevrel-phase crystals. Physical Review B     1987, 36 (4), 1952-1955. -   45. Gallenkamp, C.; Kramm, U. I.; Proppe, J.; Krewald, V.,     Calibration of computational Mössbauer spectroscopy to unravel     active sites in FeNC catalysts for the oxygen reduction reaction.     International Journal of Quantum Chemistry 2021, 121 (3), e26394. -   46. Marcus, H. L.; Schwartz, L. H., Mossbauer Spectra of FeMo     Alloys. Physical Review 1967, 162 (2), 259-262. -   47. Lu, K.; Gao, S.; Dick, R. J.; Sattar, Z.; Cheng, Y. Journal of     Materials Chemistry A 2019, 7, (11), 6038-6044. -   48. Xie, X.; Song, M.; Wang, L.; Engelhard, M. H.; Luo, L.; Miller,     A.; Zhang, Y.; Du, L.; Pan, H.; Nie, Z.; Chu, Y.; Estevez, L.; Wei,     Z.; Liu, H.; Wang, C.; Li, D.; Shao, Y. ACS Catalysis 2019, 9, (9),     8712-8718. -   49. Song, Y.; Peng, R.; Hensley, D. K.; Bonnesen, P. V.; Liang, L.;     Wu, Z.; Meyer III, H. M.; Chi, M.; Ma, C.; Sumpter, B. G.;     Rondinone, A. J. ChemistrySelect 2016, 1, (19), 6055-6061. -   50. Zhu, D.; Zhang, L.; Ruther, R. E.; Hamers, R. J. Nature     Materials 2013, 12, (9), 836-841. -   51. Watt, G. W.; Chrisp, J. D. Analytical Chemistry 1952, 24, (12),     2006-2008. -   52. Wang, M.; Liu, S.; Qian, T.; Liu, J.; Zhou, J.; Ji, H.; Xiong,     J.; Zhong, J.; Yan, C. Nat Commun 2019, 10, (1), 341. -   53. Gallenkamp, C.; Kramm, U. I.; Proppe, J.; Krewald, V.     International Journal of Quantum Chemistry 2021, 121, (3), e26394. -   54. Marcus, H. L.; Schwartz, L. H. Physical Review 1967, 162, (2),     259-262. -   55. Duan, G.; Chen, Y.; Tang, Y.; Gasem, K. A. M.; Wan, P.; Ding,     D.; Fan, M. Advances in electrocatalytic ammonia synthesis under     mild conditions. Progress in Energy and Combustion Science 2020, 81,     (100860), 1-30. -   56. Yongjun Leng, Guang Chen, Alfonso J. Mendoza, Timothy B. Tighe,     Michael A. Hickner, and Chao-Yang Wang. Solid-State Water     Electrolysis with an Alkaline Membrane. Journal of the American     Chemical Society 2012 134 (22), 9054-9057. 

1. An iron-containing Chevrel phase material, containing iron and Mo₆S₈ clusters.
 2. The iron-containing Chevrel phase material of claim 1, having a ratio, x, of iron to Mo₆S₈ cluster of 0<x≤6.
 3. The iron-containing Chevrel phase material of claim 1, having a ratio, x, of iron to Mo₆S₈ cluster of 2≤x≤4.
 4. The iron-containing Chevrel phase material of claim 1, excluding those compounds of the formula Fe_(x)Mo₆S₈ where 1≤x≤2.
 5. The iron-containing Chevrel phase material of claim 1, having a ratio, x, of iron to Mo₆S₈ cluster of x=4.
 6. The iron-containing Chevrel phase material of claim 1, having a formula Fe_(x)Mo₆S₈, wherein 2≤x≤4.
 7. The iron-containing Chevrel phase material of claim 1, having a formula Fe₄Mo₆S₈.
 8. The iron-containing Chevrel phase material of claim 1, having a formula Fe_(x)M_(y)Mo₆S₈, wherein M is at least one metal selected from the group consisting of elements of Groups 1-15 and the Lanthanide series. 0<x<6, 0≤y<6, and x+y=1 to
 6. 9. The iron-containing Chevrel phase material of claim 1, having a formula Fe_(x)M_(y)Mo₆S₈, wherein M is at least one metal selected from the group consisting of elements of Groups 1, 2, 3, 7, 10, 11, 12, and the Lanthanide series. 0<x≤4, and x+y=2 to
 4. 10. An electrode, comprising: (1) a conductive substrate, and (2) the iron-containing Chevrel phase material of claim 1, on the conductive substrate.
 11. The electrode of claim 10, wherein the iron-containing Chevrel phase material has a formula Fe₄Mo₆S₈.
 12. The electrode of claim 10, wherein the conductive substrate comprises carbon.
 13. A system for producing ammonia electrochemically from N₂ and water, comprising: (a) a working electrode, comprising (1) a conductive substrate, and (2) the iron-containing Chevrel phase material of claim 1, on the conductive substrate, (b) a counter electrode, and (c) an ion-conductive separator between the working electrode and the counter electrode.
 14. The system of claim 13, wherein container comprises: (i) a cathode compartment, with the working electrode in the cathode compartment, and (ii) an anode compartment, with the counter electrode in the anode compartment, wherein the cathode compartment and the anode compartment are separated by the ion-conductive separator.
 15. The system of claim 14, further comprising: a cathode electrolyte solution in the cathode compartment, in contact with the working electrode, and an anode electrolyte solution in the anode compartment, in contact with the counter electrode. 16-18. (canceled)
 19. A facility for manufacturing ammonia, comprising: (A) the system of claim 15, (B) a power source, (C) a water feed, and (D) an ammonia dispenser. 20-21. (canceled)
 22. A method of producing ammonia electrochemically from N₂ and water, with the system of claim 15, comprising: passing electricity through the working electrode and the counter electrode, to produce ammonia from N₂ and water.
 23. The method of claim 22, having a Faradaic efficiency of at least 12%.
 24. (canceled)
 25. The method of claim 22, having a mass-normalized yield rate of at least 50 μg/h per mg of catalyst. 26-29. (canceled)
 30. A method of producing ammonia electrochemically from N₂ and water, with the facility of claim 19, comprising: passing electricity through the working electrode and the counter electrode, to produce ammonia from N₂ and water. 